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ABSTRACT 

Context. The recent release of the First Fermi-LAT Source Catalog solidified the predominant association of extragalactic y-ray 
emitters to active galaxies, in particular blazars. A tight connection between AGN jet kinematics and y-ray properties has been argued 
for, attributing the energetic emission from active galaxies to their highly relativistic outflows. 

Aims. We investigate the Caltech-Jodrell Bank flat-spectrum (CJF) sample to study the connection between AGN jet kinematics and 
their y-ray properties. The high number of sources included in the sample, in addition to the excellent kinematic data available, allows 
us to investigate the origin of y-ray emission in AGN. 

Methods. We identify the CJF sources detected in y-rays (by Fermi-LAT and EGRET). We use y-ray luminosities and the available 
VLBl kinematic data to look for correlations between y-ray and kinematic properties, as well as for differences between AGN classes 
(quasars, BL Lacs, radio galaxies). We also check the kinematics of the TeV sources in the CJF. 

Results. 21.8% of the CJF has been detected in the y-ray regime. We find the detectability of BL Lacs significantly higher compared 
to quasars, y-detected sources show a wider apparent jet velocity distribution compared to the non-detected ones, but the maxima of 
both distributions are at similar values. No strong link between y-ray detectability and fast apparent jet speeds is found. A tentative 
correlation is found between y-ray luminosity and maximum apparent jet speeds, stronger for BL Lac and y-variable sources. We find 
non-radial jet motions to be important to y-ray emission. We suggest two-zone, spine-sheath, models as a possible explanation to our 
results. We find 2 out of 4 CJF TeV sources show superluminal jet speeds, in contrast to previous studies. 

Key words. Galaxies: statistics - Galaxies: active - Galaxies: nuclei - Galaxies: jets - Gamma rays: galaxies 



1. Introduction 

Active galactic nuclei (AGN) are known for their emission over 
a broad range of the electromagnetic spectrum, reaching up to 
energy scales of lO'^eV (e.g.. Punch et al. 1992; Neshpor et al. 
1998). This high energy end of the AGN spectrum is most prob- 
ably created by inverse Compton scattering (either of the in- 
cident synchrotron photon field, self-scattering (SSC), of pho- 
tons arising from the accretion disk, or of cosmic microwave 
background) from a population of relativistic electrons in the 
AGN jets (e.g., Blandford & Icke 1978; Ghisellini et al. 1985; 
Maraschi et al. 1992; Sikora et al. 1994). Acceleration mech- 
anisms, including acceleration from electric fields and acceler- 
ation in shocks, can boost the energy of the electron popula- 
tion enough so that y-ray emission can be efficiently produced 
through inverse Compton scattering. The variability observed at 
these short wavelengths adds a further layer of complexity (e.g., 
Hoyle 1966) and might require special geometrical and physical 
configurations to be explained (e.g., Georganopoulos et al. 2005; 
Tavecchio & GhiselHni 2008). 

Most y-ray radiation is shielded away from us by the Earth's 
atmosphere, except for very high energy (TeV) photons that 
produce atmospheric particle showers that can be detected and 
analyzed by Cherenkov telescopes. The launch of y-ray satel- 
lites like the Compton y-ray Observatory satellite (Gehrels 
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et al. 1993; carrying the Energetic Gamma Ray Experiment 
Telescope, EGRET; Kanbach et al. 1988), the AGILE satel- 
lite (Tavani et al. 2008), and the Fermi y-ray Space Telescope 
(Atwood et al. 2009) allows a detailed study of the y-ray sky. In 
particular, the Fermi telescope offers the possibility of a contin- 
uous all-sky monitoring (using the Large Area Telescope, LAT) 
at much higher resolution and sensitivity than those of previous 
instruments. There are 1451 sources included in the first source 
catalog of the Fermi-LAT (Abdo et al. 2010a), 687 of which are 
identified/associated with active galaxies. Of these, 671 Fermi- 
LAT detected sources at high galactic latitudes (\b\ > 10°) are 
included in the first catalog of active galactic nuclei (Abdo 
et al. 2010b), associated with 709 AGN (multiple associations 
for some sources), comprising 300 BL Lac objects, 296 flat- 
spectrum radio quasars, 41 AGN of other type, and 72 of un- 
known class. 

Of particular interest is the connection between the y-ray and 
radio properties of AGN and the possibly common Doppler fac- 
tor for both the radio and y-ray emitting particles. The distribu- 
tion of radio spectral indices of the Fermi-LAT detected AGN is 
consistent with a flat spectrum, although a minor tail of steeper 
spectrum sources exists (see Fig. 28 in Abdo et al. 2010b), indi- 
cating the dominance of core-dominated sources as y-ray emit- 
ters and reflecting the importance of the beaming effect for the 
production of this high-energy emission. The existence of a pop- 
ulation of steeper spectrum sources, also detected by the Fermi- 
LAT, implies a more complex link between the y and radio prop- 
erties of AGN. 
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Lister et al. (2009b) study the kinematic properties of the 
3-month Fermi-LAT detected sources of the MOJAVE sample. 
The authors find that y-detected quasars show on average faster 
jets than their non-detected counterparts, in agreement with pre- 
vious studies of EGRET-detected sources showing y-detected 
sources to have preferentially higher Doppler beaming factors 
(e.g., Jorstad et al. 2001; Kellermann et al. 2004). They also note 
that BL Lac sources show on average lower apparent speeds but 
are nevertheless preferentially detected by Fermi-LAT They at- 
tribute this behavior to the spectral shape of BL Lacs and the 
possibility that BL Lacs have a higher intrinsic y-ray to radio 
luminosity ratio than flat-spectrum QSOs, concluding that these 
results merit further investigation. 

We use the Caltech-Jodrell Bank flat-spectrum (CJF; Taylor 
et al. 1996) sample of radio-loud AGN to test the putative con- 
nection between the y-ray properties of AGN and their jet kine- 
matics, in light of the recent release of the first Fermi-LAT AGN 
source catalog (Abdo et al. 2010b) and in the context of the stud- 
ies mentioned above. The paper is organized as follows: in Sect. 
2 we describe the CJF sample, in Sect. 3 we present the available 
data, in Sect. 4 the data is analyzed and our results are presented, 
in Sect. 5 we discuss our results in the context of previous sim- 
ilar investigations, and we finally offer a summary and our con- 
clusions in Sect. 6. Throughout the paper, we assume the cos- 
mological parameters (from the first-year WMAP observations; 
Spergel et al. 2003) Hq = 71 km s'^ Mpc-\ Qm = 0.27, and 
Qa = 0.73. 

2. The CJF Sample 

The CJF sample (Taylor et al. 1996) consists of 293 sources 
selected (see Table 1) from three different samples (for details 
see, Britzen et al. 2007b). The sources span a large redshift 
range (see Fig. 1 in Britzen et al. 2008), with the furthest ob- 
ject being at a redshift z = 3.889 (1745+624; Hook et al. 1995) 
and the closest one at z = 0.0108 (1146+596; de Vaucouleurs 
1991). The average redshift of the sample is Zavg - 1.254, with 
ZBLLacavg = 0.546, ZRG,avg = 0.554, and ZQso.avg = 1-489 for BL 
Lac objects, radio galaxies, and quasars, respectively. All the ob- 
jects have been observed with the VLBA and/or the global VLBI 
network. Each source has at least 3 epochs of observations and 
has been imaged and studied kinematically (Britzen et al. 1999; 
Britzen et al. 2007b; Britzen et al. 2008). The X-ray properties 
have also been studied and correlated with their VLBI properties 
(Britzen et al. 2007a). Finally, Karouzos et al. (2010) conduct a 
multi-wavelength study of the CJF in the context of the merger- 
driven evolution of galaxies. 



Table 1. The CJF sample and its properties. 



Frequency(MHz) 


4850 


Flux lower limit @5GHz 


350mJy 


Spectral Index 


4850 > _o 


Declination 


5> 35° 


Galactic latitude 


1*1 > 10° 


# Quasars 


198 


# BL Lac 


32 


# Radio Galaxies 


52 


# Unclassified 


11 


# Total 


293 



The CJF offers an excellent tool to study the kinematics, 
given its large number of sources and the detailed investigation 
of the jet kinematics of all these sources. In light of the findings 
of Lister et al. (2009b), we are interested in investigating the 
putative connection between the y-ray properties of the AGN 
in our sample and their jet kinematics, utilizing the extensive 
database of information akeady at our disposal (Britzen et al. 
2007b; Britzen et al. 2008; Karouzos et al. 2010c, in prep.). We 
want to disentangle the different effects that might contribute in 
defining the y-ray properties of an AGN and assess their relative 
importance. In particular, we want to test whether y-detected 
AGN indeed show faster jets than their non-detected counter- 
parts. We will also investigate the jet ridge line properties of the 
y-detected sources, in relation to those of the non-detected ones 
(a detailed treatment of the CJF jet ridge lines will be presented 
in a separate paper). 



3. Data 

The CJF sample is a statistically complete sample, under the 
selection criteria described in Sect. 2. Radio data is therefore 
available for all the CJF sources, both single dish (at several 
frequencies) and interferometric. In addition, all CJF sources 
have been observed by ROSAT (see Britzen et al. 2007b for 
details). Finally, a substantial number of CJF sources has been 
recently detected in the y-ray regime by the Fermi-LAT instru- 
ment. Below we describe the radio and y-ray data available. 

3. 1 . Radio emission 

The CJF sample (Table 1) has been most extensively studied in 
the radio regime (e.g., Taylor et al. 1996; Pearson et al. 1998; 
Britzen et al. 1999; Vermeulen et al. 2003; Pollack et al. 2003; 
Lowe et al. 2007; Britzen et al. 2007b; Britzen et al. 2007a; 
Britzen et al. 2008). Britzen et al. (2008) develop a localized 
method for calculating the bending of the jet associated with in- 
dividual components. The maximum of the distribution of local 
angles is at zero degrees, although a substantial fraction shows 
some bending (0 - 40 degrees). A few sources exhibit sharp 
bends of the order of > 50 degrees (see Fig. 13 in Britzen et al. 
2008). 

Britzen et al. (2008) present an extensive analysis of the CJF 
jet kinematics, uniformly analyzing interferometric data for each 
CJF source, identifying individual components at each epoch, 
and studying their kinematic behavior. Although the CJF sample 
consists mostly of blazars, presumably highly beamed sources, 
the kinematical study of the sample shows a large number of 
sources having stationary, subluminal, or, at best, mildly superlu- 
minal outward velocities (e.g., see Fig. 15 in Britzen et al. 2008). 
Combined with a number of sources exhibiting inward mov- 
ing components (e.g., 0600+422, 1751+441, 1543+517, Britzen 
et al. 2007b), these sources do not fit into the regular paradigm of 
outward, superluminaly moving components in blazar jets. Such 
peculiar kinematic behaviors can be readily explained by geo- 
metric effects, usually assuming a helical motion pattern com- 
bined with projection effects (e.g., Zensus et al. 1995; Steffen 
et al. 1995b; Steff'en et al. 1995a; Lobanov & Roland 2005; 
Roland et al. 2008). Most of the CJF sources exhibit flux vari- 
ability at different timescales, with several showing indications 
for quasi-periodicities in their radio lightcurves (see Karouzos 
etal. 2010 for details). 
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3.2. y-ray emission 

EGRET detected 14 CJF sources and provided an upper limit 
for 50 more (Fichtel et al. 1994; Hartman et al. 1999). We note 
that, unlike for the radio and the X-rays, the y-ray study of this 
sample is not complete since EGRET only did targeted obser- 
vations of some of the CJF sources. Fermi-LAT (the renamed 
GLAST; Atwood et al. 2009) will provide a complete study of 
these sources through its all-sky survey. 61 CJF sources are in- 
cluded in the first catalog of AGN by the Fermi-LAT, after eleven 
months of observation (Abdo et al. 2010b). Three additional 
sources are included in the third EGRET catalog (Hartman et al. 
1999) but have not, as of yet, been detected by Fermi (0804-1-499, 
2346+385, and 2351+456). 

In total, we find 64 sources (21.8% of the sample) being de- 
tected in the y-ray regime, while for 40 (13.6% of the sample) 
only an upper limit is reported from EGRET. Breaking down 
the number of detections, of the 64 y-detected sources, 24 are 
classified as BL Lacs, 32 as quasars, and 5 radio galaxies. One y- 
detected source remains unclassified. We mention here that there 
are 7 CJF sources detected by Fermi-LAT that have been clas- 
sified differently by the Fermi-LAT team (see Table A.l). Three 
of these sources do not have kinematic data and therefore do not 
influence our result. For the remaining four, we retain the classi- 
fication of Britzen et al. (2007a). Six CJF sources belong to the 
group of few sources that have been detected in the TeV regime 
(Table 3). 

Of the 61 CJF sources already detected by Fermi-LAT (Table 
A.l) we note that half of them are found to be y- variable. In total, 
31 sources (10.6%) have been detected to be variable. We note 
for the following analysis we only use real detections, excluding 
upper limits information. 

It should be noted that both the inclusion of a source in the 
LAT catalog, as well as a source being flagged as variable, de- 
pends on a number of factors, the combination of which make 
the catalog incomplete. The flux of a source (both in radio and y- 
rays), its spectral shape, and its position on the sky (e.g., near the 
galactic plane) influence its possible association with an AGN. 
Moreover, as has already been noted in Abdo et al. (2010b), in 
order for a source to be classified as variable it needs to have suf- 
ficiently high y-flux. Subsequent catalogs will surely improve 
upon the current associations and consequentially increase the 
number of y-bright AGN. 



4. Analysis 

As we described in the previous sections, we are interested in 
investigating a possible link between jet kinematics and y-ray 
properties in the CJF sample. There are a couple of points that 
should be addressed before we do that. A first effect concerns 
the redshift distribution of our sources. Quasars show a redshift 
distribution peaking around z=1.3, while BL Lacs maximum is 
around z=0.2. In Fig. 1 we show the redshift distribution of the 
CJF sources detected in y-rays. The distribution follows the dis- 
tribution of the sample, with the highest y-detected source at 
redshift z = 3.044. Secondly, the CJF sample contains predomi- 
nantly quasars, with only a few percent of the sample being BL 
Lac objects and radio galaxies. This is a result of the selection 
criteria of the sample. 75% of the CJF BL Lacs have been de- 
tected by Fermi-LAT, while only 17% of the QSOs and 9.6% of 
the RGs are included in the first year AGN catalog of the tele- 
scope (Abdo et al. 2010b). Taking into account the number of 
QSOs at z < 1 (56 sources), we get a y-ray detectability of 16% 



for the CJF quasars, as expected, significantly lower than BL 
Lacs. 



Radio Galaxies 
Quasars 
BL Lacs 



N 8 



0.7 



1.4 2.1 



2.8 



Fig. 1. Redshift distribution of the Fermi-LAT detected (1st 
Catalog; Abdo et al. 2010b) CJF sources. The three differ- 
ent classes of AGN are shown: BL Lacs (dashed empty line), 
quasars (continuous empty line), and radio galaxies (continuous 
filled line). 



4.1. y-ray luminosities 

y-ray observations are based on photon counts per energy bin, 
unit surface, and unit time for each source. We are interested 
in translating the photon fluxes (at the energy interval of 1- 
100 GeV) to luminosities, in order to account for the distance 
dependence. We follow Thompson et al. (1996) to calculate y- 
ray luminosities as follows: assuming an energy range (Ei,E2), 
then the integral photon flux is described as. 



F(AE) 



r 

Je, 



AE-"dE, 



where A is a normalization constant that can be expressed in 
terms of the integral flux F and a is the spectral index. For the 
Fermi-LAT instrument, a nominal value of 1-100 GeV is used 
for the energy range probed. Then, the energy flux is calculated 
from the following formula: 



S{AE) 



1 



El-" ■ 



72-0' 



2-0- 



pi-" . 



E\-" 



where F is the photon flux, the directly observed quantity (mea- 
sured in photon counts per unit surface and time), and S is the 
energy flux (measured in ergs per unit surface and time). We use 
the median photon indices given in Abdo et al. (2010b) for each 
source. It should be noted that, within this energy range, a break 
in the spectrum might exist. This would lead to an overestima- 
tion of the luminosity, especially for the QSO sub-sample (e.g., 
Abdo et al. 2010c). 

Having calculated the sources luminosities, we can calculate 
average and median values for our sample (shown in Table 2. It 
should be noted that comparison of average and median values 
for non-Gaussian distributions is problematic. For the sake of 
continuation and uniformity with previous studies of the same 
topic we calculate both average and median values, although it 
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should be seen as only indicative and interpreted in tandem with 
results from statistical tests like the K-S or Student's test. 

Taking into account only QSOs at redshifts smaller than z= 1, 
we can calculate the average and median values (assuming a nor- 
mal distribution for the samples studied here, the uncertainty for 
the median values is calculated as 1.253 of the standard error of 
the mean) for this sub-sample (47.09 + 0.11 and47.19 + 0.14erg, 
respectively; the values are given in the log (vLy) form). It can 
be seen that QSOs are consistently more luminous than the BL 
Lacs. 

Table 2. Characteristic statistic values concerning the y-ray lu- 
minosity (in logarithmic scale, measured in erg/s) of the CJF 
sources detected in the y-ray regime by Fermi-LAT. 



log (vL,)y 


All 


QSO 


BL 


Var 


Non-Var 






All 


z< 1 










51 


31 


9 


14 


30 


21 


Average 


47.26 


47.45 


46.09 


46.32 


47.41 


46.83 


Error 


0.07 


0.09 


0.11 


0.13 


0.09 


0.09 


Median 


46.63 


46.91 


47.19 


45.44 


46.73 


46.60 


Error 


0.09 


0.11 


0.14 


0.16 


0.11 


0.11 


Max 


48.58 


48.58 


47.50 


47.15 


48.58 


47.64 


Min 


43.67 


45.33 


45.33 


43.67 


43.89 


43.67 



Notes. We use the flux between 1 and 100 GeV to calculate the y-ray 
luminosity. We give average and values with uncertainties, maximum, 
and minimum values for the whole sub-sample, QSOs, QSOs at z < 1, 
BL Lacs, variable, and non-variable sources. Given the small number 
of RGs (5 sources) detected, we do not calculate separate statistics for 
that sub-sample. 



The luminosity of a source at longer wavelengths might play 
a role in deciding its y-ray properties. To that end, we com- 
pare the luminosities of the y-detected and non-detected QSOs at 
z < 1, at 5 and 30 GHz (single dish), at optical (V band), and in 
the soft X-rays data from Britzen et al. 2007b, Taylor et al. 1996, 
and references therein). The y-detected sources are consistently 
more luminous in the radio and optical regime, but are fainter in 
the X-rays. As was described previously, up to soft X-rays the 
emission is thought to be produced by the synchrotron mech- 
anism and therefore y-detected sources show a stronger syn- 
chrotron component than their non-detected counterparts. This 
might in turn be linked to the putative inverse Compton process 
usually employed to explain the production of y-ray emission in 
AGN jets. Conversely, y-ray detection of a source implies that 
the inverse Compton hump of its SED is shifted towards higher 
energies, compared to non-detected sources. It is therefore to be 
expected that y-detected sources are actually weaker in the soft 
X-rays than their non-detected counterparts. A larger sample of 
y-detected sources, along with a closely matched (in terms of 
luminosity and redshift) control sample is required to test this 
scenario. 

Concerning the relative importance of the synchrotron and 
inverse Compton components, we calculate the y-to-radio lumi- 
nosity ratio for the y-detected CJF sources (see Fig. 2). 

In Fig. 2 we show the distribution of this ratio for quasars 
and BL Lacs (including both variable and non-variable sources). 
It can be seen that BL Lacs show marginally higher y-to-radio 
ratios than quasars, with the BL Lac distribution peaking around 
1.075 compared to the quasar one peaking around 1.055. The 
fact that we include both variable and non-variable sources 
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1.04 1.05 1.06 1.07 1.08 1.09 
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Fig. 2. Histogram comparing the distributions of y-to-radio lu- 
minosity ratios for BL Lacs and QSOs. We use y-ray luminosi- 
ties derived from the Fermi-LAT photon fluxes given in the first 
Fermi-LAT source catalog (Abdo et al. 2010a) and radio lumi- 
nosities at 5 GHz derived from single dish energy fluxes from 
Taylor etal. (1996). 



might influence our results. Given the observational bias that 
sources classified as y-variable are on average brighter than the 
non-variable ones and that there are more variable BL Lacs than 
non-variable ones, the effect observed in Fig. 2 might in part be 
due to the same observational bias discussed previously. A two 
sample K-S test for BL Lacs and quasars is inconclusive at to 
whether their y-to-radio radio distributions are drawn from dif- 
ferent parent distributions. An obvious caveat of this comparison 
is that y and radio observations are not contemporaneous. 

4.2. Apparent VLBI Jet Component Velocities and y-ray 
emission 

We characterize the kinematics of the CJF sources by the max- 
imum observed component velocity, liapp.ma.x, of each source. 
This maximum component velocity depends on both the orien- 
tation of the source (viewing angle) and the intrinsic properties 
of the jet itself. We investigate whether the distribution of maxi- 
mum apparent component speeds differs between y-detected and 
non-detected sources. For all CJF sources with available redshift, 
Britzen et al. (2008) have calculated the apparent total compo- 
nent velocities for all the identified components in their VLBI 
jets. We identify the component with the maximum apparent 
speed in each source and we plot the distribution of these max- 
imum apparent speeds in Fig. 3, for y-detected sources (black 
line) and non-detected ones (grey blocks). 

In the upper panel of Fig. 3, we plot the distribution of the 
maximum apparent velocities for the 191 CJF sources that have 
not been detected at y-ray wavelengths (grey blocks) and the 45 
sources that have (black line). The maximum of the distribution 
for the non-detections is found in the [3.5,7] bin. For the detec- 
tions the maximum lies in the [0,3.5] bin. As can be seen in the 
inset of the upper panel (inset plots are normalized to surface 
area one, as we are interested in the relative distributions of the 
two sub-samples), assuming smaller bins, the maximum of the 
y-detected sources breaks down to two maxima in the [0,2] and 
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Fig. 3. The distribution of the maximum apparent velocities fiapp,max for sources that have been detected in the y-ray regime (black 
line) and for those that have not (grey blocks). We show the distributions for all sources (upper panel), quasars (middle panel), and 
BL Lacs (lower panel). The insets show the normalized to surface area unity distributions for each case and for apparent speeds up 
to 30c. We use the kinematic data from Britzen et al. (2008). 



[4,6] bins. The non-detected sources distribution is peaked in the 
[4,6] bin. y-detected sources have a more extended distribution, 
towards higher velocities, showing a possible secondary peak in 
the [10,12] and [12,14] bins. A K-S test between y-detected and 
non-detected sources /3app,max distributions does not give a con- 
clusive answer (93.7% confidence that the two sub-samples are 
different). 

In the middle panel of Fig. 3 we show the same plots as 
before but only for those sources classified as quasars (same 
notation as before). The distribution here is markedly differ- 
ent from before. Both populations (y-detected and non-detected) 
have their distribution shifted towards higher velocities, both of 
them peaking in the [4,6] bin. The y-detected QSOs show a pos- 
sible secondary peak in the [10,12] bin (secondary to primary 
ratio of 0.71) and again appear to show a larger fraction of the to- 
tal number of sources at higher velocities. Non-detected quasars 
have a possible secondary peak in the [8,10] bin (secondary to 
primary ratio of 0.8). The highest velocities for quasars are found 



in sources that have not been detected at y-ray wavelengths. A 
K-S test gives a low probability (of 95.4%) that the two samples 
are significantly different. Similarly, in the lower panel of Fig. 3 
we show the distributions (absolute and normalized) for BL Lac 
objects. Both y-detected and non-detected BL Lacs have their 
distribution maxima in the [0,2] bin, with y-detected BL Lacs 
having a considerably more extended distribution, reaching to 
higher velocities, compared to their non-detected counterparts. 
Compared to quasars, the BL Lac velocity distribution is shifted 
towards lower values. A K-S test for the distribution of Pappjnax 
for y-detected quasars and BL Lacs gives a 97.7% confidence 
that they are drawn from different parent samples. It should be 
noted however that the number of y-detected BL Lacs with avail- 
able redshift and kinematic information is small and therefore 
our analysis for BL Lacs is probably affected by low number 
statistics. 

We also calculate average and median values of the maxi- 
mum apparent velocity for y-ray detected CJF sources and for 
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those not detected (9.0 ± 0.8 and 8.8 + 1 .0, compai-ed to 8 . 1 ± 0.4 
and 6.4 ± 0.5, respectively). A Student's t-test is inconclusive 
(49% for the null hypothesis). A further note concerns the red- 
shift distribution of the two sub-samples. The y-detected sources 
have a lower average redshift (0.967) compared to their non- 
detected counterparts (1.175). It is known that there appears 
to be a dependence between apparent velocity and redshift (as 
noted by Cohen et al. 2007; Britzen et al. 2008; Lister et al. 
2009a; Karouzos et al. 2010c, in prep.). This effect reinforces 
the difference in average and median values seen between j- 
detected and non-detected sources and implies that y-detected 
sources indeed have higher maximum apparent velocities. 

Finally, for the y-ray detected sources, we distinguish be- 
tween quasars and BL Lac objects. We find that quasars exhibit 
considerably higher average and median maximum apparent ve- 
locities (10.8 + 1.0 and 10.2 + 1.2, respectively, compared to 
6.1 ± 1.5 and 2.9 ± 1.9). Redshift effects might again be influ- 
encing our results. We select those QSOs at redshifts lower than 
1 and calculate jiapp,max average and median values. We find that 
the sub-sample of local QSOs shows, within the statistical errors, 
the same average values as the BL Lacs (6.9 ± 1.2) but consider- 
ably higher median value (6.1 ± 1.5). 

In Fig. 4 we compare the fiapp,max distribution of variable (27 
sources; average redshift z^.g = 0.94 + 0.09) and non-variable 
(19 sources; average redshift Zay^ - 0.94 + 0.10) y-detected CJF 
sources. In the following we refer to y-ray variability, not taking 
into account possible variability in other wavelength regimes. 
Non-variable y-detected sources show a more extended distri- 
bution than the variable ones, reaching the highest velocities 
(~30c). Both distributions show their main maximum in the 
[0,3.5] bin. Variable sources show a possible secondary maxi- 
mum in the [10.5,14] bin (with a secondary to primary ratio of 
0.86). A K-S test for the two distributions does not provide a 
conclusive result. 



samples show similar redshift distributions, therefore we do not 
expect any redshift effect influencing our result. 

In Fig. 5 we compare the fiapp.max distributions of y-detected, 
variable quasars and y-detected, variable BL Lacs (left panel) 
and y-detected variable quasars and y-detected non-variable 
quasars (right panel). These two sub-samples (variable sources 
and quasar sources) are the largest of the different y-detected 
sub-samples (e.g., non-variable sources, BL Lac sources, etc.) 
and therefore are chosen to check differences between quasars 
and BL Lacs, and variability, in an isolated manner 
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Fig. 5. Left panel: distributions of the maximum apparent veloc- 
ities fiapp.nmx for y-variablc BL Lacs (black line) quasars (grey 
blocks). Right panel: y-variable (black line) and non-variable 
(grey blocks) quasars. We use the kinematic data from Britzen 
et al. (2008). 
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Fig. 4. Distribution of the maximum apparent velocities I3app,max 
for y-detected sources that are variable (black line) and for 
those that are not (grey blocks). We use the kinematic data from 
Britzen et al. (2008). 



As in the previous cases, we again calculate the statistical 
properties of the two sub-samples. Variable y-detected sources 
show similar I3app,max average value to their non-variable coun- 
terparts, within the statistical errors (7.9 + 0.8 and 8.9 ± 1.7, 
respectively). When checking the median values however, the 
variable sub-sample shows substantially higher value than the 
non-variable one (8.6 + 1.0, compared to 4.6 + 2.1). Both sub- 



We find that quasars (left panel of Fig. 5) show the maximum 
of their distribution at higher values than BL Lacs ([3.5,7] bin 
compared to [0,3.5]). Given the small number for each group of 
sources, a K-S test is not applied. For variable and non-variable 
sources the comparison is not so straightforward (right panel 
of Fig. 5). Non-variable y-detected quasars show a main max- 
imum in the [0,3.5] bin, lower than their variable counterparts. 
However, they also show a more extended distribution reach- 
ing higher Pappjnax values than their variable counterparts. A K- 
S test gives inconclusive results (23% probability that the sub- 
samples are the same). Small number statistics, as well as the 
biased detection of bright sources as variable, affect our results 
and therefore we can not give a definitive answer concerning the 
relative importance of the two effects. There is evidence that y- 
ray emitting, variable quasars show statistically higher apparent 
velocities. 

The above analysis implies that there are two effects, pos- 
sible dependent, that correlate with the jet kinematic properties 
of y-emitting AGN, (1) y-variability and (2) the classification 
of the source (BL Lac or quasar). The latter can be associated 
either to the difference in un-beamed luminosity between BL 
Lacs and quasars, or possibly difference in viewing angles (e.g., 
Lahteenmaki et al. 1999; Hovatta et al. 2009). For the number of 
sources available here, we can not robustly decouple these two 
effects. Given the available data, we do not find any significant 
connection between y-ray detection and fast moving jet compo- 
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nents, as has been argued by other authors. We rather see that 
the Pappmix distribution is more strongly dependent on the type 
of object, i.e., BL Lac or quasar classification (see Sect. 5 for a 
discussion on this). 

4.3. y-ray Luminosities and Apparent Jet Component 
Velocities 

In the previous section we probed the putative connection be- 
tween the apparent brightness of y-rays and the jet kinematics 
-i.e., the beaming mechanism-, as reflected in the differences be- 
tween the j3app,nmx distribution of quasars and BL Lacs. Another 
way to approach this is by looking for a possible direct correla- 
tion between the apparent velocities measured in the AGN jets 
and their y-ray luminosity. This is shown in Fig. 6. 

As can be seen in Fig. 6, there appears to exist a correla- 
tion between the /3app,max of a source and its y-ray luminosity. 
Variable sources appear to cluster closer to the implied trend, 
with non-variable QSOs deviating the most. We find Spearman 
correlation coefficients of 0.68 and 0.77 (both at a significance 
of > 99.999%), respectively. Given the mutual dependence of 
luminosity and jet component apparent velocity to the redshift 
of a source, we also calculate the Pearson product moment par- 
tial correlation coefficients i(Ly /3„pp ,„a ^ ,z) ' . While for the whole 
sample we get a relatively low partial correlation coefficient 
(0.43 at a significance of 99.8%), the correlation for variable 
sources persists (0.68 at a significance of > 99.9%). Given the 
expected degree of scatter in the data, the correlation coefficients 
combined with the calculated significance imply that the trend 
seen in Fig. 6 is indeed true. We note that some of the appar- 
ent velocities show relatively large errors, therefore, in the right 
panels of Fig. 6, we plot only sources for which /3app,max ^ 3cr. 
We also separate variable sources from non-variable ones. It be- 
comes clear that the non-variable sources show the most scatter 
For the variable sources, we calculate a Spearman correlation 
coefficient of 0.69 (at a 99.1% significance level), lower than 
the coefficient we got when fitting all the sources. Calculating 
the partial correlation coefficient for the same sub-sample gives 
a smaller correlation coefficient (0.61 at a 98.1% significance). 
We also investigate the same correlation for individual classes of 
objects. The strongest correlations are seen for both y-variable 
BL Lacs (partial correlation coefficient of 0.68 at 95.1% sig- 
nificance) and quasars (partial correlation coefficient of 0.69 at 
99.9% significance). The differences seen between the different 
classes of AGN, as well as between variable and non-variable 
sources reveal a complicated picture. We shall discuss the ro- 
bustness and implications of our results more extensively in Sect. 
5. 

4.4. y-ray and jet ridge line properties 

The currently accepted paradigm for the jet kinematics of flat- 
spectrum sources (i.e., core-dominated AGN) includes superlu- 
minaly outward moving components, usually interpreted in the 
context of a specific projected geometry combined with relativis- 
tic effects, due to the intrinsically high speeds of the bulk flow. 
Recent detailed kinematic studies of the parsec-scale jets of BL 
Lac objects (e.g., 1803+784, Britzen et al. 2010a; 0716-H714, 
Britzen et al. 2009; 0735-H178, Britzen et al. 2010b) have how- 



' For this we use the Web tool: Wessa P., (2008), Partial 
Correlation (v 1.0.4) in Free Statistics Software (v 1.1. 23- 
r6). Office for Research Development and Education, URL 
http://www.wessa.net/rwasp_partialcorrelation.wasp/. 



ever revealed a rather different kinematic scheme for their jets: 
BL Lac jet components are predominantly stationary with re- 
spect to the core but change their position angle significantly, 
essentially reflecting an important transverse component in their 
movement. In addition, their jet ridge lines, defined as the line 
that linearly connects the projected positions of all components 
at a certain epoch, show significant temporal evolution, at times 
forming very wide flow funnels. In an accompanying paper, we 
shall present a statistical investigation of the jet ridge lines of the 
CJF sources (Karouzos et al. 2010c, in prep.). We are interested 
in investigating a possible correlation between the jet ridge line 
properties of a source and its appearance at y-ray wavelengths. 
We shall therefore briefly outline the method used to analyze the 
CJF jet ridge lines. We can define three measures that probe both 
the radial and transverse motion of the jet ridge line of a source, 
(1) the jet ridge line width, dP, (2) the jet ridge line width evo- 
lution, AP, and (3) the linear evolution of the jet ridge line, A£. 
The dP is defined as the position angle difference between the 
components with the maximum and minimum position angle at 
a given epoch, measured in degrees (see Fig. 7 for an example). 
The AP is derived between two successive epochs as the dif- 
ference between the jet ridge line widths at these epochs (mea- 
sured in degrees per unit time). The A{ is derived as the sum 
of the linear displacements, along the vector of their motion, of 
all jet components, across all available epochs (measured in par- 
sees per unit time and per component) and resembles an average 
component apparent speed across all epochs. 



aean map. Array: BFHKLMNOPS 

0700+470 «t 4.991 GHz 1995 Sep 04 dP = 9max"9min 
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Fig. 7. Example of the definition of the jet width for the source 
0700+470. With arrows we show the components with the min- 
imum and maximum position angles at that epoch. The continu- 
ous line represents the jet ridge line of 0700+470 (see the text for 
a definition) at the same epoch. Map from Britzen et al. (2007b). 



We compare these three measures of the jet ridge line prop- 
erties for y-detected and non-detected sources. We find that y- 
detected sources show significantly wider jet ridge lines, both 
in average and median values (16.9 + 1.0° and 11.1 + 1.2°, re- 
spectively), compared to the non-detected ones (13.6 + 0.4° and 
9.3 ± 0.5°). A Student's t-test gives a 97.6% probability that the 
two average values are significantly different. The difference be- 
tween the two distributions can be also clearly seen in the upper 
panel of Fig. 8: the maximum of the jet ridge line width distribu- 
tion of the y-detected (non-variable) sources is shifted to higher 
values (in the [5,10] bin), compared to their non-detected coun- 
terparts (which have their maximum in the [0,5] bin). 

Furthermore, y-detected sources are found to show stronger 
evolution of their widths, both in average and median values 
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Fig. 6. Maximum apparent component velocities for all y-ray detected CJF sources as a function of their y-ray luminosity, for all 
sources (left panel) and for high significance values of Papp.max (right panels; see text for details). We differentiate between y-variable 
(filled symbols) and non-variable sources (open symbols), and between the different AGN classes, i.e., BL Lacs (circles), quasars 
(squares), and radio galaxies (diamonds). We use the kinematic data from Britzen et al. (2008). 
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Fig. 8. Histograms of jet ridge line width, dP, and width evo- 
lution, Af , for y-detected variables sources (continuous black 
line), y-detected non-variable sources (dashed black line), and 
y-ray non-detected sources (grey blocks). Given the large dif- 
ference in absolute numbers, we plot histograms normalized to 
surface area unity. We use data from Karouzos et al. (2010c, in 
prep.). 



(4.7 ± 0.5 deg/yr and 2.5 + 0.6 deg/yr, respectively), com- 
pared to their non-detected counterparts (3.13 + 0.18 deg/yr and 
2.26 ± 0.22 deg/yr). A Student's t-test confirms at a 98.1% con- 
fidence that the two averages are significantly different. In the 
lower panel of Fig. 8 we show the normalized distribution of the 
width evolution values for y-detected and non-detected sources. 
The y-detected distribution is peaked in the [2,3] bin, compared 
to the non-detected one, peaking in the [0,1] bin. 

We can also distinguish between variable and non-variable 
y-detected CJF sources and investigate their jet ridge line width 
evolution in this context. Variable sources are found to have sig- 
nificantly stronger evolving jet ridge line widths, at a 4cr signif- 
icance level, both in average and median values (6.7 +1.0 and 
4.4 ± 1.2 deg/yr, compared to 2.4 ± 0.3 and 2.2 ± 0.4 deg/yr). 
A Student's t-test confirms this result (at a 98.8% confidence). 
In Fig. 8 (lower panel) we plot the distribution of y-detected 
variable and non-variable sources. The non-variable sources dis- 
tribution is fairly confined to lower values, whereas variable 
sources extend up to the highest values of width evolution. This 
implies that viewing angle changes at parsec scales of AGN jets 
are linked to variability in the y-ray regime, reflecting a possibly 
evolving jet at even smaller scales (connected to the timescales 
associated with y-ray variability). 

We finally compare how the jet ridge line evolves, in linear 
terms, in y-detected and non-detected sources, essentially this 
time looking for a possible link between y-ray brightness of a 
source and its average jet component apparent speed, instead of 
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the maximum apparent speed. As we already discussed previ- 
ously, the redshift distributions of the two sub-samples are fairly 
similar. We do not expect therefore any redshift-induced effects, 
related to linear distances, to influence our results. We find that 
both in average and median values, the two sub-samples show 
similar values of linear evolution (0.47 + 0.03 and 0.40 ± 0.04 
pc/yr/comp, respectively, for the y-detected sources, compared 
to 0.440 ± 0.020 and 0.377 ± 0.025 pc/yr/comp for the non- 
detected ones). 

The same behavior is seen when comparing between y- 
detected variable and non-variable sources. It should be noted 
however, that the median value for y-detected variable sources is 
fairly higher than the one for non-variable sources (0.43 + 0.04 
pc/yr/comp compared to 0.291+0.025 pc/yr/comp, respectively). 
This discrepancy between average and median values implies a 
large scattering within our data that might be posing some limi- 
tations to the robustness of our results. Aside from that, we con- 
firm our finding from before, i.e., that the link between fast ap- 
parent jet speeds and the y-ray detectability of a source is ques- 
tionable. It appears that, although there are certainly indications 
that higher speed sources are preferentially y-ray emitters, some 
other effect might play a more important role in defining the 
y-ray properties of a source. Our results concerning the width 
and width evolution comparison between y-detected and non- 
detected sources imply that non-radial motion -i.e., motion trans- 
verse to the jet symmetry axis- is important with regard to y-ray 
production. This is probably tightly connected to the viewing an- 
gle effect, also presumably reflected in the kinematic differences 
between different AGN types. 

4.5. TeV sources 

28 extragalactic sources have been identified to be emitters in the 
TeV regime (for an updated list of these sources and references, 
see http://tevcat.uchicago.edu/). Six CJF sources (5 BL Lac ob- 
jects, 1 radio galaxy) are observed and detected in TeV. All are 
variable in y-rays and optical, five in the infrared, while five of 
them are variable in radio and four in the X-rays. 

Extremely high energetic photons coming from AGN (TeV 
sources, see Table A.l) can be a result of several different mech- 
anisms, some of them taking place in the AGN jets (for a dis- 
cussion on high energy y-ray emission from AGN see e.g., von 
Montigny et al. 1995). We check the kinematics of the TeV 
sources in the CJF (Table 3; data from Britzen et al. 2008). We 
give the misalignment angle between the parsec and the kilo- 
parsec scale (as calculated by Britzen et al. 2007a). We also 
list the bending of the parsec scale jet (as calculated by Britzen 
et al. 2008) for all identified components. Finally we give the 
maximum apparent total jet component speed, /3app,,nax- We find 
no apparent correlation between the bending of the jet and the 
y-ray flux. 2 out of 5 sources exhibit superluminal motion in 
their VLBI jets. Previous studies of TeV emitting high-frequency 
peaked BL Lacs exhibit slow jet speeds (e.g.. Finer et al. 2008; 
Britzen et al. 2009). 0219-H428 and 2200-H420 do not belong to 
that category of objects. 

5. Discussion 

It is interesting to compare our results with similar studies using 
different samples of AGN. As we already mentioned in Sect. 1, 
a number of authors have claimed a close connection between 
y-ray bright AGN and large [iapp,max (e.g., Jorstad et al. 2001; 
Kellermann et al. 2004; Lister et al. 2009b). Lister et al. (2009b), 
in particular, examine the kinematic properties of the MOJAVE 



sources (Lister et al. 2009a) in light of the, then recently, pub- 
lished 3-month bright AGN Fermi-LAT list. In short, they find 
that (1) the y-detected MOJAVE quasars have the peak of their 
Papp,max distribution around 10- 15c, quite a higher value com- 
pared to the non-detected ones (~5c). They also find that (2) y- 
variable sources show higher apparent speeds. Finally, although 
BL Lacs in their sample have lower redshifts and slow median 
jet speeds (6c), (3) they are preferentially detected by Fermi- 
LAT. Similarly to Lister et al. (2009b), Jorstad et al. (2001) and 
Kellermann et al. (2004) also conclude that y-detected AGN 
have significantly higher apparent speeds than their non-detected 
counterparts. 

The results presented here are somewhat different. Although 
we do find a possible secondary maximum around 10- 15c for the 
y-detected sources, the primary maximum of both distributions 
appears to be around 5c, with y-detected sources actually show- 
ing a double maximum in the [0,2] and [4,6] bins. The fiapp^max 
distribution of the sources is more strongly dependent on the 
classification of the source as a quasar or a BL Lac, and whether 
the source appears variable, rather than it being detected at y- 
ray wavelengths. Given the larger number of sources in the CJF, 
unlike Lister et al. (2009a), we can treat the different AGN pop- 
ulation separately. In particular, as we showed in the previous 
section, quasars show higher apparent jet speeds than BL Lacs, 
regardless of their y-ray properties. A common feature for y- 
detected sources is that their [iapp,max distribution appears more 
extended, with larger percentage of the sources at higher speeds, 
compared to those of non-detected sources. Coming to point (2) 
-i.e., that y-variable sources show higher apparent jet speeds- 
our results agree, as we showed for the case of variable and 
non-variable y-detected quasars. However this does not come 
as a surprise in light of our recent findings from the statistical 
analysis of the CJF jet ridge lines. In short, we find that vari- 
able sources (throughout the electromagnetic spectrum) show 
stronger linear evolution of their jet ridge lines, a measure essen- 
tially reflecting their apparent jet speeds (Karouzos et al. 2010, 
in prep.). Therefore, the higher apparent speeds of y-detected 
variable sources can probably be linked to the variability, rather 
than to the y-ray brightness itself. 

Finally, concerning point (3), we also find that BL Lac ob- 
jects are preferentially detected by Fermi-LAT, even after cor- 
recting for the different redshift distribution of quasars and BL 
Lacs. Lister et al. (2009b) argue that BL Lacs might be preferen- 
tially detected by Fermi-LAT because of their flatter spectrum, 
i.e., their higher y-to-radio ratio, compared to quasars. As was 
shown in Sect. 4.1, this is indeed true. 

A further point of interest, that has not so far been inves- 
tigated, concerns the apparent correlation seen between the y- 
ray luminosity of a source and its maximum apparent jet speed. 
As we showed in Fig. 6 such a correlation exists, with vari- 
able sources falling closer to the implied trend. We need to con- 
template whether the inclusion of more sources with potentially 
lower y-ray luminosities would destroy this correlation. Looking 
at the CJF sources not yet detected by Fermi-LAT, 35 of them 
show significant (> 3cr) apparent speeds above 10c, while only 
10 show speeds larger than 20c. Of these 10 sources only 3 are 
classified as high quality component fits (see Britzen et al. 2008 
for details). Together with the high median redshift of the high- 
speed y non-detected CJF sources, the above suggest that the 
correlation seen in Fig. 6 is robust. 

Although our results are generally in agreement with previ- 
ous studies, there are some differences that merit further inter- 
pretation. Of particular interest is the fact that we do not find 
a strong connection between fast moving components and y- 
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Table 3. Jet morphology and kinematics of TeV detected CJF sources. 



Source 


Type 


z 


APA 

(deg) 


Bend 

(deg) 


Plot,max 


y-rays 

(10"**^/? cm-- 5"') 


0219+428 


BL 


0.444 


8 


24 

4 

1 


14 ±5 


2.49 ±0.10 


0316+413 


G 


0.018 


20 


52 
85 
83 
24 


0.98 ± 0.20 


1.73 ± 0.08 


0716+714 


BL 




75 


3 
7 
3 




1.31 ± 0.07 


1101+384 


BL 


0.031 


13 


2 


0.19 ± 0.07 


2.61 ± 0.10 


1652+398 


BL 


0.034 


83 


10 
40 
10 
3 


0.90 


0.83 ± 0.06 


2200+420 


BL 


0.069 


30 


25 
16 
47 


3.2 ±0.3 


0.71 ±0.06 



Notes. Columns (l)-(3) give the lAU name, the type, and the redshift of the source. Col. (4) gives the misalignment angle between the pc and 
kpc-scale jet (from Britzen et al. 2007a), Col. (5) gives the bending of the identified components of the pc-scale jet (for details see Britzen et al. 
2008), Col. (6) gives the maximum apparent total speed for the VLBl jet of the source, and Col. (7) gives the source's y-ray flux between 1 and 
100 GeV (Tabel A.l). 



detected sources, as opposed to previous works. As we already 
discussed in Sect. 1, one of the models proposed for the origin of 
-y-ray emission in AGN involves a spine-sheath geometry, where 
a high velocity spine is embedded in a slower moving sheath. 
According to these models, the ultra-relativistic spine gives rise 
to the y-ray emission that we observe. Assuming that the jet as 
a whole becomes gradually more opaque at lower observing fre- 
quencies, we could then argue that, given the lower frequency 
that the CJF sample has been observed in, we are only prob- 
ing the (transversely) outer regions of the jet and therefore we 
only recover the slower motions linked to the sheath of the jet. 
Conversely, the MOJAVE sample (15 GHz) as well as the sam- 
ples used by both Jorstad et al. (2001) (42, 22, and occasionally 
15 and 8.4 GHz) and Kellermann et al. (2004) (15 GHz) are all 
at higher observing frequencies and thus more sensitive to the 
inner layers of the jet. It is therefore plausible that these studies 
have recovered the faster speeds linked to the ultra-relativistic 
spine and the production of the observed y-ray emission. An ad- 
ditional elfect related to the higher frequency observation of the 
MOJAVE programme is the distances from the core probed. At 
15 GHz, the higher resolution of the VLBI observations allow 
the investigation of smaller spatial scales, closer to the core, and 
therefore possibly more closely connected to y-ray production 
processes. Finally, given the limited number of epochs available 
for some of the CJF sources, combined with the low observing 
frequency, it can be argued that the fastest moving components 
might have been missed or blended together with slower compo- 
nents in the CJF kinematic study, therefore resulting to, on av- 
erage, slower component speeds for the whole sample. Detailed 
studies of individual sources at both higher temporal and spa- 
tial resolutions (e.g., 1803+784, Britzen et al. 2010a; 0716+714, 
Britzen et al. 2009; etc.) have shown that this is not the case (if 
not showing the presence of the reverse elfect). 

An alternative scenario would be that indeed y-ray emis- 
sion, although boosted in high-Doppler jets, is not necessarily 
coupled to fast jet components speeds. Given that Lister et al. 



(2009b) use the three month bright source list for their study, 
it is possible that only the y-brightest, hence the most strongly 
beamed, sources are detected. Therefore there is a strong bias to- 
wards sources with higher Doppler factor jets. Including sources 
from the first, 11 -month, Fermi-LAT catalog allows us to probe 
lower-flux sources (given the higher sensitivity achieved after 
11 -months of observations) and therefore the bias towards the 
most highly beamed sources is lifted. This scenario however is 
unlikely since (1) by using only the sources included in the 3- 
month bright sources Fermi-LAT list and our CJF apparent speed 
data, we still recover a low-velocity component in the Papp,ma.x 
distribution of the y-detected sources and (2) Savolainen et al. 
(2010), after calculating Doppler and Lorentz factors of Fermi- 
LAT detected AGN (3-month list), find a notable absence of 
sources at the smallest (and therefore most highly beamed) co- 
moving viewing angles. Distinguishing between either scenarios 
will require a larger sample of y-detected AGN with uniform, 
multi-frequency kinematic data available. 

A last possibility that might be affecting the kinematic prop- 
erties of the y-detected sources pertains to the relativistic ef- 
fect of sources observed at a viewing angle inside the 1 /F cone, 
where F is the Lorentz factor of the flow, having lower appar- 
ent speeds. Hovatta et al. (2009) calculate median values for 
Lorentz factors and viewing angles for different AGN samples. 
The authors find for BL Lac objects a median Lorentz factor of 
F = 10.29 and a median viewing angle of = 5.24, while for flat- 
spectrum quasars they find 16.24 and 3.37, respectively. These 
values suggest that both sub-classes of AGN are found inside the 
1 /F cone. Given a certain combination of intrinsic speeds and 
viewing angles, sources could occur that are beamed enough to 
be picked up by Fermi-LAT but at the same time exhibit slower 
apparent jet speeds. The fact that the same population of sources 
is not seen in the MOJAVE sample makes this explanation de- 
batable. 

A further layer of complexity is introduced when we con- 
sider the y-ray variability of most of the y-detected sources. If 
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the y-ray emission from AGN is of transient nature, related to 
a relativistic shock propagating down the pc-scale jet, this im- 
plies that many sources are missed from such an investigation, 
simply because they are observed during their quiescent state. In 
addition to this, the non-contemporaneity of the kinematic and 
y-ray data introduces further noise to plots like Figs. 3 and 6. At 
present we can not sidestep these two caveats. 

The last point that needs to be discussed concerns the con- 
nection between the jet ridge line width to the y-ray properties 
of a source. Pushkarev et al. (2009) study the jet opening an- 
gles of the MOJAVE sample again in the context of the 3-month 
bright source list of the Fermi-LAT. The authors find that y- 
bright sources have larger apparent opening angles, but in their 
co-moving frame both y-bright and faint sources show similar 
opening angle distributions. This is in turn interpreted as evi- 
dence that y-detected sources are seen at smaller viewing angles 
than the non-detected ones. Our results agree with these findings, 
as we find that y-detected sources exhibit wider jet ridge lines 
compared to their non-detected counterparts. It should however 
be noted that preferential detection of BL Lacs at y-ray wave- 
lengths combined with our results from the statistical analysis of 
the CJF jet ridge lines (i.e., that BL Lacs have apparently wider 
jet ridge lines; Karouzos et al. 2010c, in prep.) might introduce 
a spurious effect in such an analysis. It is also interesting to note 
that, according to our findings from Sect. 4, y-detected sources 
also exhibit large changes in their jet ridge line widths, larger 
than the non-detected ones. That would introduce a noise factor 
when comparing the width distributions of y-detected and non- 
detected sources. This in turn implies that y-detected sources 
might in reality be even wider than inferred by the distribution 
shown in Fig. 8. 

It is interesting to contemplate on how a two-zone model 
might accommodate for the apparently important non-radial 
component of the component trajectories, as implied by the 
significantly larger jet ridge line width evolution shown by j- 
detected sources. A helical structure of AGN jets combined with 
the smaller viewing angles expected for y-bright sources might 
offer a plausible explanation. Such a geometry has been used to 
explain non-ballistic trajectories for a number of AGN jets (e.g., 
3C345, Steffen et al. 1995b; 1803+784, Steffen et al. 1995a; 
1633H-382, Liu et al. 2010). One could speculate that a heli- 
cally structured sheath, enclosing an ultra-relativistic (straight) 
spine flow, could give rise to the kinematic behavior observed in 
the CJF sources, while explaining the y-ray properties of these 
sources. Given the still small number of available CJF sources 
detected by Fermi-LAT, it is, for the time being, not possible to 
decouple and assess the individual importance of these effects. 

6. Conclusions 

In the previous sections we have investigated the connection of 
the y-ray properties of the CJF sources to the morphologic and 
kinematic properties of their jets. In summary we find: 

- 21.8% of the CJF sample is detected at y-ray wavelengths 
(either from EGRET or the Fermi-LAT; three EGRET asso- 
ciations not included in the 1 1 month catalog) 

- BL Lacs appear to be preferentially detected in the y-ray 
regime. Taking into account the difference of the redshift dis- 
tributions of QSOs and BL Lacs in our sample, we still get 
fairly different detection ratios between the two classes (16% 
and 75%, respectively). 

- after calculating the y-ray luminosities of both QSOs and BL 
Lacs, and taking into account the redshift effects of our flux- 



limited sample, we find that QSOs appear more luminous at 
y-ray wavelengths than BL Lacs. 

- y-detected sources (regardless of classification) show the 
peak of their Pappjnax distribution at similar values as their 
non-detected counterparts but show a more extended distri- 
bution towards higher jet component speeds. A K-S test re- 
sults in an inconclusive answer, giving a 93.7% confidence 
that the two sub-samples are significantly different. 

- when considering QSOs and BL Lacs separately, we still find 
that y-detected and non-detected sources show fairly simi- 
lar distributions, with y-detected sources showing more ex- 
tended distributions towards higher values. 

- comparing between y-detected QSOs and BL Lacs, after ac- 
counting for redshift effects, both sub-samples show on av- 
erage and within statistical errors the same apparent speeds. 

- we find a tentative correlation between fiapp.max and y-ray lu- 
minosity. The correlation is stronger for BL Lac objects and 
for y-variable sources, with non-variable QSOs deviating the 
most from the implied trend. 

- we find that y-detected sources have significantly wider jet 
ridge lines than their non-detected counterparts. We also find 
that y-detected sources show stronger jet ridge line width 
evolution than non-detected ones. 

- we find no significant difference in terms of hnear evolution 
of the jet ridge lines between y-detected and non-detected 
sources. 

- we find no direct link between highly bent jets and Te V emis- 
sion. Furthermore, we note that 2 out of 4 TeV CJF sources 
with fiapp,max information show superluminal speeds, unlike 
previous studies. 

From our analysis it becomes clear that there is a number 
of factors influencing whether a source is luminous in the y-ray 
regime or not. Although it is tempting to think in terms of y-loud 
and quiet objects, the picture is surely more complicated. As has 
been recently demonstrated for the radio divide (radio-loud and 
radio-quiet objects) there is a number of sources that actually 
seem to populate an intermediate space, implying a rather con- 
tinuous distribution. It is possible that with more sensitive y-ray 
telescopes a greater number of y-ray emitting sources will be re- 
covered. Of course the question still remains as to what it is that 
differentiates the y-quiet (or -faint) sources from those already 
detected by the Fermi-LAT and EGRET missions. Our analysis, 
combined with previous studies on this question, indicates that 
the viewing angle and y- variability are what make some sources 
to be y-ray luminous and others not. However, as has been shown 
by Savolainen et al. (2010), the picture may be more complicated 
than that. Given the relatively weak link between jet apparent 
speeds and y-detected sources that we find in this work, as well 
as the different picture arising at higher observing frequencies, 
a spine-sheath configuration scenario could offer a plausible ex- 
planation, where the most energetic emission - coming from an 
ultra-relativistic flow- originates in the spine of the AGN outflow 
and is therefore partly or fully obscured at lower observing fre- 
quencies. Finally, the link between the width and width evolution 
of the jet ridge lines of our sources to their y-ray properties im- 
plies that transverse, non-radial, motions in the jet might be im- 
portant in this context, probably also related to the y-variability 
detected for almost half of the y-detected AGN in our sample. 
It is obvious that a larger sample of y-detected AGN with com- 
plete, multi-frequency kinematic data will allow us to investi- 
gate the above effects in a more robust manner. Furthermore, 
the investigation of both individual objects and statistical sam- 
ples (like the CJF and the MOJAVE, e.g., Chang et al. 2010; 
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Chang et al. 2010, in prep.) with detailed modehng of their SEDs 
will surely shed light to the actual processes producing the y-ray 
emission. 
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Appendix A: y-ray properties of the CJF sample 



Table A.l. CJF sources detected in the y-ray regime by Fermi-LAT. 
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z 


T? 
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1 -2 -1\ 

phcm ^ s ) 
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log(vLv)y 6\o%{vL^)y 
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0.6 


0.3 
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0.13 
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0.22 


N 


4.899 
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Q 
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0.7 


0.3 
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0.18 


45.33 


0.19 


N 
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0.646 
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Q 
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0.6 
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0.03 


47.33 


0.03 
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Q 
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1 
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0.13 
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0.17 


Y 
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0218H-357 


G(Q) 


0.936 


6.4 


0.5 


2.33 


0.04 


47.25 


0.03 


Y 


n/a 


n/a 


0219-1-428 


BL 


0.444 
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1 
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0.09 


Y 
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1.07 
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Y 
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17.3 
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Y 
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2 
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Y 
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BL 
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1 
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2.03 


0.18 
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N 
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Q 
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0.6 


0.3 


2.73 


0.17 
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0.22 


N 
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Q 
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0.5 


2.32 


0.04 


47.23 


0.04 


Y 


n/a 


n/a 
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0.3 
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0.14 


N 
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BL 
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0.7 
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Y 
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BL 
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N 
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Y 
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BL 
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0.03 


Y 
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Q 
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Y 
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Q 


2.18 
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N 
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Q 
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14 
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48.58 


0.02 


Y 
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Q 
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0.3 
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N 
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BL 
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0.12 


N 


n/a 


n/a 
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Q 
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0.02 


N 
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n/a 
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BL 
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0.26 


N 
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Q 
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0.2 
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0.15 
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0.17 


N 
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1.59 
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Q 
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N 
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Q 


1.12 


1.1 


0.3 


2.48 
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N 
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Y 
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0.08 
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N 
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BL 
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Q 
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Y 
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Q 
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0.14 


N 


n/a 
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BL 
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N 
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BL 
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2.45 


0.15 
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N 
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BL 
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2.52 


0.18 
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N 
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BL 


0.069 


7.1 
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Y 
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Notes. Columns (l)-(3) give the lAU name, type, and redshift of the source. In brackets is given the alternative classification trom the Fermi-LA'l' 
team, where applicable. Column(4) gives the y-ray flux (detections are from the first Fermi-LAT catalog, Abdo et al. 2010a), Col. (5) gives the 
uncertainty for the y-ray fluxes. Cols. (6) and (7) give the y photon index and its uncertainty. Cols. (8) and (9) give the calculated y-ray luminosity 
(in the log vLy form and in logarithmic scale) and the respective uncertainty. Col. (10) denotes sources with detected variability in y-rays, and Cols. 
(11) and (12) show the maximum apparent jet speed for that source and its uncertainty. 
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